Boson stars: Chemical potential and quark condensates 
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We study equilibrium solutions of a boson star in mean field approximation using a general 
relativistic framework. Dynamics of the constituent matter in such a star is described by a scalar 
field. The statistical aspect of the matter equilibrium is investigated by incorporating the effect of 
finite chemical potential in equation of state. We analyze the two possible situations where the scalar 
field is either a composite one -suitable for describing diquark-condensates or a fundamental. We 
derive a generalized set of Tolman-Oppenheimer-Volkov (TOV) equations to incorporate the metric 
dependence of chemical potential in general relativity. It is demonstrated that the introduction of 
finite chemical potential can lead to a new class of the solutions where the maximum mass and 
radius of a star change in a significant way. In general the density profile has node like structures 
due to introduction of finite chemical potentials. It is shown that these nodes can be avoided by 
applying constraints on the values of the central pressure and central chemical potential. This in 
turn can reduce the parameter space available for the stable solutions. We also discuss the case 
when the boson star is made of diquark condensates. It is shown that when the self-interaction of 
the condensate is negligible, the typical energy density of the star is too high to have a boson star 
of diquark condensates. Our results indicate, even after the effect of self-interaction is incorporated, 
the stable equilibrium can occur only in a low-density and high-pressure regime. 

PACS numbers: 



I. INTRODUCTION 



The boson stars are very fascinating objects as their 
self-gravity is not balanced by the degeneracy pressure 
like the other compact stars such as a white dwarf or 
a neutron star, but the Hciscnberg uncertainty play a 
crucial role in their stability (see Ref. [l[ for a recent 
review). Wheeler was first to consider the boson star in 
1955 when he studied the stable equilibrium generated by 
the self-gravitating photons called geons Q • Later Kaup 
analyzed a new class of boson stars as a self-consistent 
solutions of Einstein-Klein-Gordon equations in Ref. Q . 
Since then the most of the boson star studies focus on 
the the stars made of self-gravitating scalar-field [l|, 3] ■ 
This is because there are good reasons to believe that 
there exists a some fundamental scalar field in nature 
even though no experiment has detected it so far. The so- 
lutions found by Kaup in Ref. [1] were shown to describe 
the ground state of boson stars without using any perfect- 
fluid approximation or some approximation to equation 
of state by Ruffini and Bonazzola Q. However in this 
work the self-interaction of the scalar-field was ignored. 
Subsequently it was shown by Colpi et. al. Q that in- 
clusion of even a very small interaction term can signifi- 
cantly alter the solutions obtained in Refs. [1,[H. It turns 
out that the inclusion of a small interaction term in the 
scalar field Lagrangian give rise to an additional pressure 
that would help the Heisenberg uncertainty to counter 
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the self-gravity. This results in a boson star with much 
larger mass than ones reported in Refs. Also the 

issues of dynamical stability of boson stars have been rig- 
orously analyzed by applying perturbative methods [7[. 
Also are studies that show that there can be a boson star 
in the galactic center Q. At present a very intriguing 
situation persists in the field: on the one hand there is 
no observational evidence which can indicate that a bo- 
son star exists, on the other hand their existence seem to 
be a reasonable consequence of well-tested physical the- 
ories like general relativity and Klein- Gordon equations. 
This should be a good enough reason, in our opinion, to 
further investigate the properties of boson stars. 

Boson stars with a repulsive self-interaction mediated 
by vector mesons within the effective theory framework 
were studied only recently It is found that mass- 
radius curves satisfy the scaling properties for arbitrary 
boson masses and interaction strengths. There exist an 
another novel possibility of having a boson star with a 
composite scalar field. Such a composite scalar field can 
possibly arise when the superconducting colour quark 
matter undergo a transition from the BCS phase to Bose- 
Einstein condensation phase via a so called BCS-BEC 
transition [13, [HI- It should be noted that BCS regime 
of the colour superconducting phase is well studied using 
the techniques of perturbative QCD. But in BCS-BEC 
crossover can occur in low density and low temperature 
regime, where the coupling constant is large and the per- 
turbative techniques are not available. One may use the 
methods of effective theory to describe the matter in the 
strongly coupled regime. In this picture, instead of a 
repulsive vector-meson exchange one may have repulsive 
self- interaction of the scalar field describing the (diquark) 
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condensates [T3| • A diquark-BEC phase may occur at low 
density regime somewhere between colour neutral nuclear 
matter and the free quark phases. This requires an ad- 
ditional constraint namely - the size of the condensates 
or the coherent length ^ should not be larger than the 
inter-particle spacing i.e. n~^l'^ , where n is the number 
density. 



In this work we analyze the equilibrium solutions of a 
boson star in a low temperature finite chemical poten- 
tial regime. As far as we know boson stars with a finite 
chemical potential are not properly considered in the lit- 
erature. Majority of study do not focus on the statistical 
aspect of the self-gravitating boson gas but they rather 
appeal to find a specific equilibrium (quantum) solution 
of the Einstcin-Klcin-Gordon system. In the context of 
quark star we would like to mention that possibility of di- 
quark star has already been considered in the literature 
However the effect of Bose- Einstein condensation 
was not considered in this work. This effect can be sig- 
nificant in a low temperature finite chemical potential 
regime [lo| . Moreover, in general relativity temperature 
and chemical potential arc local functions of space and 
time and therefore they depend upon the metric. The 
metric dependence of these quantities can be specified 
by using the well-known Tolman conditions for a ther- 
mal equilibrium in the gravitational field [l^ • Again 
this aspect was not considered in the literature on boson 
or diquark stars. In this work we derive a generalized set 
of Tolman- Oppenheimer-Volkov equations by incorporat- 
ing the effect of metric dependence of chemical potential. 
This has resulted in an additional differential equation 
describing the space-time evolution of the chemical po- 
tential coupled with the TOV [l5| equations. 



In a very high density system such as a boson star, the 
introduction of finite chemical potential can rise a ques- 
tion. However in the standard model, for example, quark 
numbers are conserved and one can introduce chemical 
potential for them. In addition for a BCS-BEC crossover 
kind of transitions chemical potential play a very impor- 
tant role Chemical potentials are known to alter 
the global minimum of a scalar-field. In this work have 
demonstrated that when the scalar field mass m is com- 
parable to the chemical potential /i, the mass-radius rela- 
tions and the maximum mass of the star are significantly 
changed. 



In what follows we first derive the generalized set of 
TOV equations and study their numerical solutions by 
providing the equation of state as an input. The equation 
of state is derived in the mean-field approximation. We 
use for the parameter values of the diqua rk p hase in a 
manner consistent with that given in Ref. M. 



II. FORMALISM 



Lagrangian density for the scalar-field can be written 



as 



with 



- (1) 

where, the potential is defined as 

The energy- momentum tensor r^jy($) written from this 
as follows: 

V = [a^$ta,$ + a^<i>a,$t] (3) 

The equations of motion of the scalar field is given by 



where 
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We are considering the 'standard' form of a spherically 
symmetric and static space-time metric [T6| given by 



-B{r)dt^ + A{r)dr^ + r^{d0^ + sin^ 



'■), (6) 



where, we have used the units in which speed of light 
c = 1. Next to study the equilibrium configuration, we 
use the stationarity ansatz ^{r,t) = ^(r)e~*'^*. Now the 
energy density p and pressure p can be obtained from the 
energy-momentum tensor: 
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And finally the equation of motion for the scalar field 
gives. 
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It must be noted here that while writing equations ([?]■ 
[5]) we have dropped the spatial derivative terms of the 
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scalar field 
A^^. Here A 



It can be shown that is of the order of 



is an extremely large number 



(> 10'^°) for the parameter regime of the interest. This 
is equivalent to the mean field approximation (9^$ = 0). 
Under this condition one can use the TOV equation to 
describe a self-gravitating bose star @. 

Using equations ^HQ^ one can rewrite p and p 
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where, po — p — and p — p/ po is the dimension- 
less pressure. Equation (fTTI) is equivalent to the EOS con- 
sidered in Refs. P, 0, Il7( in the limit ^ = 0. In the limit 
of vanishing self- interaction i.e g — >■ 0, we have p = 0. 
In this limit self-gravity of the star will be balance d by 
the effect of Heisenberg uncertainty principle only [3lll8|. 
For high values of the p the last term in parenthesis of 
equation (11) dominates over the other terms, the EOS 
approach the relativistic limit i.e. p = 3p. 

The form of functions A{r) and B{r) in equation ([6]) 
can be determined from the Einstein's field equations. 



Rfj.L' — -Qn.ijR — —SttGT, 
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The equations ([T5)) and are known as the TOV 
equations [lH|, which has to be solved together with the 
equations of state. But in our case p is a function of chem- 
ical potential /x, which in turn depends on the metric. 
Metric dependence of p can be understood using equa- 
tion (jl7p and thermodynamic Gibbs-Duhem relation. 
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and can be written as [ij 
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Finally the we get the desired relation with the help of 
equation as. 
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So generalized stellar structure equations in presence of 
chemical potential are given by equations (jl6p . (|18p and 
(PT|). These equations are to be solved simultaneously 
with equation ([TT]), to find out the properties of stellar 
configurations. 



With the metric defined previously we get the relevant 
equations from (IT^ as, 
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where ' denotes differentiation with respect to r. With 
A{r) = (1 — 2GM{r)/r)~^ , we get the stellar structure 
equations as. 
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Hydrostatic equilibrium expressed as [li 
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can be used to rewrite the equation (fT5|) as 
G[p + p] [M + Anr^p] 
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III. RESULTS AND DISCUSSION 

Before solving stellar structure equations numerically 
let us consider the following qualitative arguments. First 
consider the case when there is no self-interaction term 
and also no chemical-potential i.e. g = and p = 0. In 
this case the critical mass and radius of the star given 
^ and M = M^;^„^j,/m Next, if the boson mass 
is around lOOGeU, the corresponding critical mass of the 
star would be M = ~ 10^% ~ lO'^^AfAr 

and its radius R ^ IO'^^Rn where, Mat and Rn are 
mass and radius of a neutron star. Consequently, the 
density of such a boson star would be 10^^ times that of 
neutron star. Thus for a diquark-condensate with mass 
(0.5 — l)GeU [l^, the density of the corresponding star 
would be 10*"'^ times of that of a neutron star. This is 
an extremely large value of the matter density for a di- 
quark star to exist. However, things can change if the 
self-interaction term is nonzero. It is known that for the 
case g (and p — 0), the scalar-field | $ | scales as 
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This condition is 

to satisfy even for a very small values of g. Comparing 
the quadratic and quartic terms in the potential U , we 



get 



A 



1*1^ 



-pla 



0(1). The energy density p will 



be about m'^M^i^^^^./A. Thus the energy density would 
corresponds to that of a non-interacting boson if mass is 

rescaled as m — > to-i/4-. From this the maximum mass 
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can be found [1, [l^ to be 
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This is also matches with the estimate critical mass ob- 
tained numerically in Ref. @. In our case, in the flat 
space, the classical potential has non-trivial minimum 
at I $ p= ^ ~" . We have assumed that > 
and for /i^ > we have condensation. In this case 
the mass scaling will be different and it is given by 

m m/VA - l) and the maximum mass may be 

given by VA ^^"^'""'' - 1^ This suggest that the 

maximum mass should decrease if the chemical poten- 
tial is decreasing. Since the chemical potential term here 
depend on space-time the above is only a rough estimate 
of the maximum mass. It also must be noted that the 
above arguments are not vahd if < 1. 

The validity condition ^ < n^^/'^ for the condensate 
description. can be found from the gap equation 

M 



A = 2fiexp — 



(23) 



for n « 400 — 500MeV and the corresponding running 
coupling constant 3.56. However, at this juncture the 
above equation should not be taken too seriously as it 
is based upon the arguments of the perturbation theory. 
The BCS-BEC transition is supposed to take place in the 
high coupling regime. The quark-condensation condition 
in the flat-space is However, in the case 

of a boson star it needs to be locally satisfied. 

Next, we write the generalized TOV equations in a 
dimensionless form. We have already introduced the di- 
mensionless variables for the pressure p, density p and 
chemical potential /i. We calculate the mass M in the 
dimensionless unit M = M/Mq where Mo is the solar 
mass. Radius of the star (R) is also computed in terms 
of ^ = R/Ro, with Ro = GMo = I AT km. Radius of 
the star can be defined as the distance from the center 
where p becomes zero. We solve equations p6l [TSl [2T|) 
by specifying various values of the central pressure Pc 
and the chemical potential Pc, while M at the center is 
considered to be zero. 

From the form of equation (11) it is clear that, 
the energy density could be negative if the values Pc 
and Pc are chosen arbitrarily. Therefore it is required 
to constrain their values by imposing the condition p > 0. 

/5 > condition and two branches of solutions : 
The condition (on the initial conditions) implies that 
there are following two branches of the solutions: 
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FIG. 1: The graph shows the solutions for initial conditions 
of pressure pc and chemical potential pc for pc > 0. 
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The solid curve in figure (1) represents the first branch 
given by equation (25), which is valid the values of chemi- 
cal potential defined in the range 

1 < A < The curve 
with the dotted line represents the branch given by equa- 
tion (26). From figure (1) it is clear that solutions of the 

equation (25) is valid even in the region 1 < p < 

Figure (2) shows the density profile as function of f for 
various values of the central pressure and the chemical 
potential p = 1.5. The solution indicates the nodes in 
the profiles. However, this is a unphysical behavior and 
it can be avoided by imposing the no-node condition on p. 

No-Node condition : We need to make sure that 
p'(r = 0) < inside the star. Analytically this condi- 
tion can be represented as, 



^Vm) U 1 + Vm - m) p'm (26) 



4(1 + y/m) p{0)p'{0)<0 



The upper most curve in figure (2) is plotted by in- 
corporating the no-node condition on the initial values 
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FIG. 2: For a given chemical potential the density profile can 
have nodes. The upper most curve is obtained by imposing 
the no-node condition. 



FIG. 3: Maximum Mass of the star versus initial pressure for 
different values of flc is plotted. The peak of the graph in 
each case denotes the highest maximum mass configuration. 



of pressure and chemical potential. Thus the implemen- 
tation of the no-node condition shows that for a given 
chemical potential there exist a lower bound on the cen- 
tral pressure. If one specify the central pressure below 
this bound then the nodes develops in the density pro- 
files. Our numerical study indicates that if the initial 
values of the pressure and chemical potential arc chosen 
just above the two branches of 'p > condition' as shown 
in figure (1), the density profile develops the nodes. One 
can avoid the nodes for higher values of the central pres- 
sure Pc, but this will also increase the central density as 
implied by equation (fTTj). It should be noted that the 
nodes in our analysis arise solely due to the introduction 
of the chemical potential and therefore their origin is dif- 
ferent than the nodes in the wave functions reported in 
earlier literature of boson-star (for example @). 

Figure (3) shows graph of the star-mass as a function 
of the central pressure for different values of the chemi- 
cal potential. The lower most curve represents the case 
of zero chemical-potential. The mass is increasing with 
increasing the central pressure (density) and it reaches a 
plateau. This behavior is in agreement with the earlier 
work. Also the value of maximum-mass is in agreement 
with Ref. @. The introduction of a finite chemical- 
potential increases the overall value of the maximum- 
mass and the plateau. This figure shows that when the 
values of fXc, increases, the maximum-mass values occurs 
at higher values of Pc- The lower values of pc in this case 
(/ic > 1) are ruled out due to the no- node condition. If 
this condition is violated then the maximum-mass con- 
figurations would occur at much lower pc- And in this 
case the maximum-mass would be much higher than the 
/Ltc = 1 case. Thus the consistent solutions with finite 
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FIG. 4: Maximum Mass of the star versus Radius for different 
values of /ic is plotted. The peak of the graph in each case 
denotes the highest maximum mass configuration. 



chemical-potential requires to have higher values of the 
central-pressure in comparison with the /ij, = case. 

In figure (4) we have plotted mass of the star ver- 
sus radius for different values of the chemical-potential. 
The radius is determined by the value of radial coordi- 
nate r (= i?), where the density vanishes. The metric- 
dependence of the chemical-potential gives a new class 
of solutions. Figure (4a) shows the case of zero central 
chemical-potential, in this case there is no variation in 
the chemical-potential as implied by equation (PT|). In 
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this case, the maximum mass value occurs at the peak 
of the curve and its value matches well with the known 
result . All the points on the left-hand side of the peak 
are unstable, as they have higher gravitational potential 
energy than the points on the right-hand side. Figure 
(4b) represents /ic = 1 case, one can see that the value 
of maximum-mass has increased as compared to /.ic = 
case and more strikingly the value of the maximum-mass 
does not occur on the peak of the curve. The spiral 
like structure in the lower part of the curve represents 
two branches of then solution. In this region for any 
given R there exist two points with the same mass. The 
points on the lower-arm of the spiral have minimum en- 
ergy, while points on the upper-arm may be unstable. 
Figure (4c) shows the case with fic = 1.2. Again the 
maximum-mass does not occur on the peak of the curve 
and its value is smaller than /ic = 1 case. Radius of 
the star at the maximum-mass point is much smaller 
than the case shown in figure (4b) and thus the den- 
sity of the star is higher than = 1 case. Figure (4d) 
shows that for /ic = 1-5, the maximum-mass is smaller 
than /ic = 1.2 case and it occurs at a smaller radius. 
In this case the maximum-mass point lie in the unsta- 
ble region. The lower-arm of the spiral represents the 
stable branch of the solution. The maximum-mass point 
in the stable branch is comparable to the case shown in 
figure (4a). This points occurs when (^^)~^ = is sat- 
isfied. It should be noted that the qualitative arguments 
for maximum- mass of a star, given below equation (22), 
may not be applicable for the cases shown in figure (4b- 
4d). As these arguments do not account for the variation 
in the chemical-potential and no-node condition imposed 
on the p. 

Figures (5-6) show pressure and density as functions 
of f for the case when the central-pressure is kept fixed 
at pc = 35 and vary flc- Figure (5) indicates that the 
pressure decreases less slowly with r by increasing the 
value of chemical-potential. Thus the star radius may 
increase if the chemical-potential is increasing. The den- 
sity vs. distance plot also shows the similar behavior. In 
this case the value of central-density pc is not fixed, but 
it decreases as the chemical-potential at center increases. 
This behavior is expected from the form of the EOS ([TT|). 

Figure (7) shows the typical variation in chemical- 
potential and the metric-function B with the radius. The 
metric-function approaches the unity at the boundary, 
while the value of the chemical-potential decreases from 
its central value with the radius. However, the value of p 
remains non-zero at the boundary. The dynamics of the 
chemical-potential and B are complementary as can be 
anticipated from definition of p. 

Let us consider the case when the scalar-field repre- 
sents the diquark-condcnsates which may arise in a low- 
density and strong-coupling regime. The condition for 
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FIG. 5: The variation of the normalized pressure with 
the normalized radius for the different values of the central 
chemical-potential is shown here. The central pressure in all 
the cases is kept fixed. 
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FIG. 6: The variation of the normalized density with nor- 
malized radius for the different values of the central-chemical 
potentials is shown here. The central pressure in all the cases 
is kept fixed. 



the condensate formation is {p/m)'^/B > 1 or p^ > 1. 
One can notice from figure (7) that this condition is vi- 
olated somewhere inside the star. For the condition to 
be valid inside the entire region of the star one has to 
increase the value of pc- From the discussion above, in- 
creasing Pc would also requires to increase the central- 
pressure Pc- This may in turn increases the density as 
shown in figure (2). Clearly the energy density can not 



7 



2.0 



1.5 









— 'f 
... B 


\ / 

\ / 
\ / 

v 








/ 

/ 

/ 



2 4 6 8 10 12 14 



FIG. 7: The typical variation of the normaUzed chemical- 
potential and the metric-function B with the radius is shown 
for the given values of central pressure and chemical-potential. 

increase to a very high values otherwise the condition of 
BEC formation may not be fulfiUed. For example one 
may take the value of the chemical potential flc = 1-5 
which is consistent with [l^. According to the figure (2), 
the physical density profile should be around 30 times 
larger then the nuclear density which is not consistent 
with the BCS-BEC crossover assumption. Moreover, the 
mass-radius diagrams shown in figure (4) suggest that as 
fl increases the available parameter-space for the stable 



configuration reduces significantly. We would like to add 
here that when one considers the case of fundamental 
scalar-field, the constraint of the low energy-density is 
not required. In this situation the boson-star equilibria 
are defined over a large range of the parameters. 



IV. SUMMARY 

Thus we have studied the properties of boson-stars in 
the mean-field approximation. We derived the general- 
ized set of TOV equations to incorporate the effect of 
finite chemical potential. It is shown that chemical po- 
tential can introduce unphysical features like nodes in the 
equilibrium profile of the density which can be removed 
by imposing the no-node condition. No-node condition 
can put constraints on the physical values of the chemi- 
cal potential and the pressure at the center of the star. 
The new class of solutions of the star that we have found 
show that available parameter space for the stable so- 
lutions gets greatly reduced by increasing the chemical 
potential. Further when our analysis when applied to 
the diquark condensate stars show that for a stable solu- 
tion the density of the star needs to be too high for the 
diquark condensates to exist. 
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